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We propose a numerical algorithm for finding optimal mea- 
surements for quantum-state discrimination. The theory of 
the semidefinite programming provides a simple check of the 
optimality of the numerically obtained results. With the help 
of our algorithm we calculate the minimum attainable error 
rate of a device discriminating among three particularly cho- 
sen non-symmetric qubit states. 

Nonorthogonality of quantum states is one of the ba- 
sic features of quantum mechanics. Its deep conse- 
quences are reflected in all quantum protocols. For in- 
stance, it is well known that perfect decisions between 
two nonorthogonal states cannot be made. This has im- 
portant implications for the information processing at the 
microscopic level since it sets a limit on the amount of 
information that can be encoded into a quantum system. 
Although perfect decisions between nonorthogonal quan- 
tum states are impossible, it is of importance to study 
measurement schemes performing this task in the opti- 
mum, though imperfect, way. 

Two conceptually different models of decision tasks 
have been studied. The first one is based on the min- 
imization of the Bayesian cost function, which is nothing 
but a generalized error rate In the special case of 
linearly independent pure states, the second model - 
unambiguous discrimination of quantum states — makes 
an interesting alternative. The latter scheme combines 
the error-less discrimination with a certain fraction of 
inconclusive results 

Ambiguous as well as unambiguous discrimination 
schemes have been intensively studied over the past few 
years. In consequence, the optimal measurements distin- 
guishing between pair, trine, tetrad states, and linearly 
independent symmetric states are now well understood 
iflroJsHl^ll- Many of the theoretically discovered opti- 
mal devices have already been realized experimentally, 
mainly with polarized light Jl^-[2l[|. As an example of 
the practical importance of the optimal decision schemes 
let us mention their use for the eavesdropping on quan- 
tum cryptosystems p2| , p3[ . 

The purpose of this paper is to develop universal 
method for optimizing ambiguous discrimination be- 
tween generic quantum states. 

Assume that Alice sets up M different sources of quan- 
tum systems living in p-dimensional Hilbert space. The 
complete quantum-mechanical description of each source 



is provided by its density matrix. Alice chooses one of 
the sources at random using a chance device and sends 
the generated quantum system to Bob. Bob is also given 
M numbers specifying probabilities that i-th source 
is selected by the chance device. Bob is then required 
to tell which of the M sources {pi} generated the quan- 
tum system he had obtained from Alice. In doing this he 
should make as few mistakes as possible. 

It is well known JjJ] that each Bob's strategy can be 
described in terms of an M-component probability op- 
erator measure (POM) {IF,-}, < IF < 1, Ej n j = l - 
Each POM element corresponds to one output channel 
of Bob's discriminating apparatus. The probability that 
Bob points his finger at the A:-th source while the true 
source is j is given by the trace rule: P(k\j) — TrpjUk- 
Taking the prior information into account, the average 
probability of Bob's success in repeated experiments is 



P, 



M 



(1) 



Since the objective is to keep Bob's error rate as low as 
possible we should maximize this number over the set of 
all M-component POMs. In compact form the problem 
reads: 



maximize P s subject to constraints 
IF (I. j = 1 . . . M, 



(2) 



Unfortunately, attacking this problem by analytical 
means has chance to succeed only in the simplest cases 
(M = 2) Q, or cases with symmetric or linearly indepen- 
dent states |^|i^,[l5|,|l7],Q . In most situations one must 
resort to numerical methods. In the following we will 
use the calculus of variations to derive a simple itera- 
tive algorithm that provides a convenient way of deal- 
ing with the problem (^). This approach has already 
found its use in the optimization of teleportation proto- 
cols [p5| and maximum-likelihood estimation of quantum 
measurements [^6| . We are going to seek the global max- 
imum of the success functional P s subject to the con- 
straints given in Eq. (j^). To take care of the first con- 
straint we will decompose the POM elements as follows 
IF, = AjAj, j = 1...M. The other constraint (com- 
pleteness) can be incorporated into our model using the 
method of uncertain Fagrange multipliers. Putting all 
things together, the functional to be maximized becomes 
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C = X ^TrteAjAj} - Tr{A ]T A/Aj}, 



(3) 
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where A is a Hermitian Lagrange operator. This ex- 
pression is now to be varied with respect to M inde- 
pendent variables Aj to yield a necessary condition for 
the extremal point in the form of a set of M extremal 
equations for the unknown POM elements: £jpjHj = 
ALTj, j — 1 . . . M. For our purposes it is advantageous to 
bring these equations to an explicitly positive semidefi- 
nite form, 

II, ;;A l Pj U jPj X '. j = l...M. (4) 

Lagrange operator A is obtained by summing Eq. (Q) over 
j, 



1/2 



Pi 



(5) 



The iterative algorithm comprised of the M+ 1 equations 
(f|) and (||) is the main formal result of this paper. One 
usually starts from some "unbiased" trial POM {n^}. 
After plugging it in Eq. (||) the first guess of the Lagrange 
operator A is obtained. This operator is, in turn, used 
in Eq. (||) to get the first correction to the initial-guess 
strategy {n°} p7| . The procedure gets repeated, until, 
eventually, a stationary point is attained. Notice that 
both the positivity and completeness of the initial POM 
are preserved in the course of iterating. 

Since equations (^) and (||) represent only a necessary 
condition for the extreme, one should always check the 
optimality of the stationary point. In the following we 
will make use of the theory of the semidefinite program- 
ming (SDP) |2^] to derive a criterion of the optimality 
of the iteratively obtained POM, which turns out to be 
the well-known Helstrom condition jjj. SDP theory also 
provides alternative means of solving the problem (J2j) nu- 
merically. 

Recently, it has been pointed out j2Sj ] that many prob- 
lems of the quantum-information processing can be for- 
mulated as SDP problems. For instance, let us compare 
our original problem (^) to the SDP dual problem that 
is defined as follows: 



maximize 
Z > 0, 

TiFjZ = c;, i 



TtFqZ, 



(6) 
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where data are m + 1 Hermitian matrices Fi and a com- 
plex vector c 6 C m , and Z is a Hermitian variable. As 
can be easily checked, our problem (|^) reduces to a dual 
SDP problem upon the following substitutions: 



*b = -0&ft, /■ 011,. 



Fi = @Ti, Ci = TrTi, i=l. 

3=1 



■P 2 . 



(7) 



Here operators {Ti,i — l...p 2 } comprise an orthonor- 
mal operator basis in the p 2 -dimensional space of Hermi- 
tian operators acting in the Hilbert space of our problem: 
TiTjTfc = 5jk, j, k = 1 . . .p 2 . For simplicity, let us take 
Ti proportional to the unity operator, then all Cj apart 
from ci vanish. 

An important point is that there exists a primal prob- 
lem associated with the dual one, 



minimize c x, 

F(x) =F + J2 t XiFi > 0. 



(8) 



Here data Fi and are the same as in Eq. (^|), and vector 
x is now the variable. 

The advantage of the SDP formulation of the quantum- 
state discrimination problem is that there are strong nu- 
merical tools designed for solving SDP problems. In par- 
ticular, these methods are guaranteed to converge to the 
real solution. This might become important when the it- 
erative algorithm derived above encounters convergence 
problems. 

A primal (dual) SDP problem is called "strictly fea- 
sible" if there exists x (Z) satisfying the constraints in 
Eq. (||) [Eq. (||)] with sharp inequalities. One can easily 
check that both the primal and dual problems associ- 
ated with the quantum-state discrimination problem are 
strictly feasible. Hence we can use a powerful result of 
SDP theory saying that in this case, x is optimal if and 
only if x is primal feasible and there is a dual feasible Z 
such that 



ZF{x) 



0. 



(9) 



This condition is called the complementary slackness con- 
dition. Now, taking Xi to be the coordinates of the La- 
grange operator A in Pj basis, Xi — TrAr^, i — 1 . . . p 2 , the 
complementary slackness condition is seen to be equiva- 
lent to the extremal equation. Since the dual feasibility 
of the iteratively obtained POM elements is guaranteed 
by construction, our extremal equation becomes a neces- 
sary and sufficient condition on the maximum of P s once 
the positive semidcfinitcness of F(x) is verified. In terms 
of states and POM elements this latter condition reads: 



A - i,/< ; > 0, j = 1 . . . M 



(10) 



One perceives that complementary slackness condition 
(^) together with the criterion of optimality (|l0|) are 
nothing else than the well-known Helstrom equations 0] 
for POM maximizing the success probability ([j]). 

Let us illustrate the utility of our algorithm on a sim- 
ple, albeit nontrivial example of discriminating between 
three non-symmetric coplanar qubit states. The geome- 
try of this problem is shown in Fig. |l|. 
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FIG. 1. A cut through the Bloch sphere showing the states 
to be discriminated. 

and ^2 taken to be equal-prior states, £i = £2 = 
£/2, symmetrically placed around the z axis; the third 
state lies in the direction of x or y. A similar configu- 
ration (with ^3 lying along z) has recently been inves- 
tigated by Andersson et al. ]l7| ]. Exploiting the mirror 
symmetry of their problem the authors derived analytic 
expressions for POMs minimizing the average error rate. 
For a given angle tp the optimum POM turned out to have 
two or three nonzero elements depending on the amount 
of the prior information £. 

Our problem is a bit more complicated one due to the 
lack of the mirror symmetry. Let us see whether the 
transition from the mirror-symmetric configuration to a 
non-symmetric one has some influence on the qualitative 
behavior of the optimal POMs. Minimal error rates cal- 
culated using the proposed iterative procedure [Eqs. (^) 
and (|^)lfor the fixed angle of p ~ tt/16 are summarized 
in Fig. |. 
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FIG. 2. Average error rate (1 — P s ) in dependence on Bob's 
prior information £; ip = n/16. Regions I, II, and III are 
regions where the optimum discriminating device has two, 
three, and two output channels, respectively. 

The conclusions that can be drawn from the numerical 
results partly coincide with that of Ref. (l7). For large 
£ (region III) the optimum strategy consist in the op- 
timal discrimination between states and ^2- When 
£ becomes smaller than a certain ^-dependent thresh- 



old (region II), state '5 3 can no longer be ignored and 
the optimum POM has three nonzero elements. Simple 
calculation yields 



I + sin if cos ip 

for the threshold value of the prior. However, when £ 
becomes still smaller (region I), the optimum POM will 
eventually become a two-element POM again - the op- 
timal strategy now being the optimal discrimination be- 
tween states and ^3. This last regime is absent in the 
mirror-symmetric case. The transition between regions I 
and II is governed by a much more complicated expres- 
sion than Eq. (]TT|) , and will not be given here. We will 
close the example noting that already a few iterations are 
enough to determine the optimum discriminating device 
to the precision the elements of the realistic experimental 
setup can be controlled with in the laboratory. 

In this paper we derived a simple iterative algorithm 
for finding optimal devices for quantum-state discrim- 
ination. Utility of our procedure was illustrated on a 
non-trivial example of discriminating between three non- 
symmetric states. From the mathematical point of view, 
the problem of quantum-state discrimination is a prob- 
lem of the semidefinite programming. Such correspon- 
dence is a good news since there exist robust numerical 
tools designed to deal with SDP problems. These can 
substitute our iterative algorithm in the very few excep- 
tional cases (if there are any) where our procedure might 
suffer from the convergence problems. 
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